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We present a Next-toLeading order perturbative QCD analysis of 
world data on the spin dependent structure functions g^, gi, and gf, 
including the new experimental information on the dependence 
of y". Careful attention is paid to the experimental and theoretical 
uncertainties. The data constrain the first moments of the polarized 
valence quark distributions, but only qualitatively constrain the 
polarized sea quark and gluon distributions. The NLO results are 
used to determine the dependence of the ratio gi / Fi and evolve 
the experimental data to a constant = 5 GeV^. We determine 
the first moments of the polarized structure functions of the proton 
and neutron and find agreement with the Bjorken sum rule. 
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1 Introduction 



Next-to-Leading Order (NLO) perturbative QCD (pQCD) analyses of unpo- 
larized lepton-nucleon deep inelastic scattering (DIS) [1-3] have resulted in 
the decomposition of the structure functions into valence quarks, sea quarks 
(of each flavor), and gluons. The data upon which these analyses are based 
include the scattering of charged leptons, neutrinos, and antineutrinos off a va- 
riety of targets, including both protons and deuterons, over a large kinematic 
range in both Bjorken x and momentum transfer Q^. 

Presently data are also available for polarized DIS [4-14] . Values for the polar- 
ized structure functions g\ {x) have been measured for protons, neutrons, and 
deuterons over a reasonable region of x and with good precision. Analyses 
of the first moments of the structure functions, Fi = / gi(x)dx, have indi- 
cated that relatively little of the spin of the nucleon is carried by the quarks, 
suggesting that perhaps the sea quarks and gluons are polarized. Hence it is 
desirable to decompose the spin-dependent structure functions into contribu- 
tions from valence quarks, antiquarks, and gluons just as has been done for 
the spin-averaged structure functions. 

On the theoretical side, a full calculation of the Next-to-Leading Order (NLO) 
spin-dependent anomalous dimensions has been recently completed [15]. This 
provides for a perturbative QCD analysis of polarized DIS as a tool for decom- 
posing the structure functions [16-19]. However, the lack of polarized neutrino 
data and the limited kinematic coverage in x and of the polarized DIS data 
limits the conclusions that can be drawn. 

We have recently reported on a precision measurement of the neutron spin- 
dependent structure function at an average four-momentum transfer squared 
= 5 GeV^ in SLAG experiment E154 [13]. The two independent spectrom- 
eters used in E154 provided a possibility of studying the dependence of the 
structure function g^, and extended the kinematic range of the measurement 
beyond that of the previous SLAG experiments [7-10] to 0.014 < x < 0.7 and 
1 GeV^ < < 17 GeV^. The El 54 results presented in this Letter supple- 
ment our previously published data [13]. They currently constitute the most 
precise determination of g^". 

Of special interest for our data is the observation that the absolute value 
of g^ increases rapidly as x becomes small for x < 0.1, approximately as 
^-0.9 j-j^gj rpj^jg jg striking contrast with the assumption of Regge behavior, 
which suggests that g^ is constant or decreases in magnitude with decreasing 
X [20]. Moreover, if the observed x-dependence of g^ persists to x — 0, the 
first moment becomes unrealistically large. 

We will show that by using NLO pQGD and reasonable assumptions about the 
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relation of the polarized and unpolarized distributions, we can obtain excellent 
fits to our data which can be used to determine the first moments and F". 
Based on these fits, we evaluate what we know about the polarization of gluons 
and sea quarks. Careful attention is paid to the theoretical and experimental 
errors involved in the analysis. 



2 Formalism 



In the QCD-improved quark-parton model (QPM), the polarized structure 

function gi{x) of the nucleon is related to the polarized quark, antiquark, and 
gluon distributions Aq{x), Aq{x), and AG{x) via the factorization theorem 
[21] 



Nf 



Cq (g) {Aq + Aq) + —Cg ® AG 



with the convolution (g) defined as 



(C0q) {x,Q') = J'^c(^-^,as)q{z,Q'). 



(2) 



The sum is over all active quark flavors Nf. 

The first moments of the structure functions of the proton and neutron, Ff 
and F", allow one to test the fundamental Bjorken sum rule [22] and deter- 
mine the helicity content of the proton. The information on the x and 
dependence gives insight into the perturbative and non-perturbative dynam- 
ics of quarks and gluons inside the nucleon. Coefficient functions Cq^G(a;, a^) 
correspond to the hard scattering photon-quark and photon-gluon cross sec- 
tions and are referred to as Wilson coefficients. They are calculated in pQCD 
as an expansion in powers of the strong coupling constant as- In leading or- 
der, C^^^ = (5(1 — x) and Cq^ = according to the simple partonic picture. 

The polarized NLO coefficient functions (7^^-* and Cq^ in the modified mini- 
mal subtraction (MS) renormalization and factorization schemes are given in 
Ref. [15]. Throughout this paper, we use the fixed-flavor scheme [2,16] and set 
Nf — 3 in Eq. (1). The heavy quark contributions are included in the run- 
ning of the strong coupling constant as{Q'^) calculated to two loops [23]. For 
consistency with the evolution of the unpolarized distributions, we adopt the 
values of asiQ'^) and current quark masses from Ref. [2] that correspond to 
as{M^) = 0.109 or as{5 GeV^) = 0.237. We include the uncertainty associ- 
ated with the value of as as will be discussed below. The parton distributions 
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in Eq. (1) are those of the proton. The neutron structTirc ftmction is ob- 
tained by the isospin interchange u d, and the dcutcron structure function 
is defined as gf = {l/2){gi + gi){l — I.Bujd), where the i^-state probabihty 
ujd = 0.05 ±0.01 [24]. 

The evolution of the parton densities is governed by the DGLAP equations 
[25] 



where the index NS stands for the the non-singlet quark distributions: valence 

Auv{x,Q'^) = Au ~ Au, Adv{x,Q'^) = Ad - Ad, and the SU(3)flavor non- 
singlet combinations Aq-^^x, Q"^) = {Au + A?;,) — {Ad + Ad) and Aqs{x, Q"^) = 
{Au + Au) + {Ad + Ad) - 2{As + As). The SU(3)flavor singlet distribution 
is AJ:{x,Q^) = {Au + Au) + {Ad + Ad) + {As + As). The index 77 = 1 
refers to the evolution of the valence (charge-conjugation odd) distributions 
uy and dy, and rj = —1 refers to the evolution of the charge-conjugation 
even combinations Ag3, Agg, and AS. The splitting functions P^g and Pij are 
calculated perturbatively with the leading order functions given in Ref. [25], 
and the next-to-leading order expressions recently obtained in Ref. [15]. In 
leading order, the evolution of both types of non-singlet distributions is the 
same: P^s^~~^ = Pt^s^~^^ = P^ and the differences only appear in ncxt- 
to-leading order. Starting with a parametrization of the parton densities at 
low initial scale = 0.34 GcV^, the distributions at any value of > Qq 
are obtained using the solutions of the NLO DGLAP equations in the Mellin 
n-moment space [26,27]. The structure functions evolved in MeUin space are 
inverted back to Bjorken x space using the prescription of Ref. [27]. 

One of the conventions relevant to the interpretation of the deep inelastic 
scattering data at next-to-leading order arises form the relative freedom in 
defining the hard scattering cross sections C^^q and the singlet quark density 
AE in Eq. (1), known as the factorization scheme dependence [26,28,29]. In 
the unpolarized case, the factorization scheme is fixed by specifying the renor- 
malization procedure for the hard scattering cross sections Cg^a [26,29]. In the 
polarized case, the situation is further complicated by the freedom in the defi- 
nition of the 75 matrix and the Levi-Civita tensor in n 7^ 4 dimensions [28,30] 
in dimensional regularization [31]. The NLO splitting functions and Wilson 
coefficients are given in Ref. [15] in the MS scheme with the definition of the 75 
matrix following Ref. [31]. The specific feature of this scheme is that the first 
moment of the gluon coefficient function vanishes C^^^n = 1) = 0, and the 
gluon density does not contribute to the integral of gi. Several authors [32- 
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34] have advocated using a different factorization scheme in which the axial 
anomaly contribution — (a5((5^)/47r) 5]]^ e^AG is included in the integral Fi. 
The suggestion generated a vivid theoretical debate [28,32-35]. Such a scheme 
was referred to in Ref. [17] as the Adler-Bardeen (AB) scheme. In the AB 
scheme the total quark hehcity is redefined compared to the MS scheme 



ZTT 

AE^ = Ago(g'), (4) 

where Aqq is the proton matrix element of the SU(3)flavor singlet axial current. 
An attractive feature of the AB scheme is that ASab is independent of even 
beyond the leading order. One could also resurrect the naive QPM expectation 
AE fa 0.6 — 0.7 and explain the violation of the EUis-Jaffe sum rule if the 
product as{Q^)/^G{Q'^) turned out to be large [32-34]. 

The product as{Q'^)AG{Q'^) is scale-independent in the leading order since its 
anomalous dimension expansion starts at order q;| [36]. This implies that as 
as decreases logarithmically with Q^, AG grows as l/as{Q'^). This growth is 
compensated by the increasing (with opposite sign) orbital angular momentum 
contribution (L^) [37,38] in order to satisfy the sum rule 

1aE + AG+(L,) = 1. (5) 



The gauge-invariant and scheme-independent formulation of this sum rule has 
recently been presented in Ref. [39]. 

Another consequence is that the ambiguity in the definition of the total quark 
helicity in Eq. (4) does not vanish at infinite Q^. However, as long as the fac- 
torization and renormalization schemes are used consistently, NLO predictions 
can be made for the spin dependent structure functions and other hadronic 
processes involving spin degrees of freedom once the parton distributions are 
determined in one scheme and at one scale. 

A transformation from the MS scheme of t'Hooft and Veltman [31] to the AB 
scheme was constructed in Ref. [17]. This scheme is a simple modification of 
MS since it preserves the low and high x behavior of the coefficient functions 
and anomalous dimensions, and thus the asymptotic behavior of parton distri- 
butions is not modified. In order to demonstrate the effects of the factorization 
scheme dependence, we perform our calculations in both MS and AB schemes. 
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3 Fits 



Following Ref. [16], we make our central ansatz of parametrizing the polarized 
parton distribution at the low initial scale Qq — 0.34 GeV^ as follows: 

^f{x,Ql)=AfX^f{l-xpf{x,Ql) , (6) 



where A/ = Auy, Ady, AQ, AG are the polarized valence, sea, and gluon 
distributions (see below for the definition of AQ). and /(,x, Qq) are the unpo- 
larized parton distributions from Ref. [2]. The parametrization assumes the 
power-like asymptotic behavior of the polarized distributions at low x and 
low Q^, namely Af ~ x'^f , x — > 0, where 7/ is the sum of the polarized 
power af and the low x power of the unpolarized distribution. Since inclu- 
sive deep inelastic scattering does not provide sufficient information about the 
flavor separation of the polarized sea, we assume an "isospin-symmetric" sea 
Au — Ad = I (Au + Ad) . Under this assumption, the sea quark contribution 
to the polarized structure functions of the proton and neutron is the same: 

gf - g'l '''' = {5/9)Cg ® \l/2{Au + Ad) + l/5Asl . (7) 



Inclusive DIS does not probe the light and strange sea independently. The 
only sensitivity to the difference between Au, Ad, and As comes from the dif- 
ference in the evolution of the two types of non-singlet distributions {rj = ±1 
in Eq. (3)). Such a difference is beyond the reach of present-day experiments. 
Hence, we will parametrize a particular combination of the sea quark distri- 
butions that appears in Eq. (7): 

AQ = 1/2 ( Am + Ad) + 1 /5As . (8) 



Furthermore, we assume the x dependence of the polarized strange and light 
sea to be the same, and fix the normalization of the strange sea by 

Au + Ad K . ^ 

A. = A.- = A,^^ = ^^A«, (9) 



with the SU(3)flavor symmetry breaking parameter Xg varying between 1 and 
(where the latter choice corresponds to an unpolarized strange sea). 

The positivity constraint, |A/(a;)| < f{x), satisfied (within uncertainties) at 
the initial scale Qq, holds at all scales Q^ > Qq; it leads to the constraints 
CKf > and Pf > 0. In addition, we assume the helicity retention properties of 
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the parton distributions [40] that require^ f3f = 0. Unhke most NLO analy- 
ses [16-18], we do not assume SU(3)flavor symmetry and do not fix the normal- 
ization of the non-singlet distributions by the axial charges Aq^ = F + D and 
Agg = 3F — D, where F and D are the antisymmetric and symmetric SU(3) 
couphng constants of hyperon beta decays [41]. Thus, we are able to test the 
Bjorken sum rule. In addition, the structure functions are not sensitive to the 
corrections beyond NLO in the data range. 

The remaining eight coefficients are determined by fitting the available data 
on the spin dependent structure functions g^'^'"^ of the proton [6,8,10,11], neu- 
tron [7,13,14], and deuteron [9,10,12] with > 1 GeV^. We use either the 
results for gi or determine the structure functions at the experimental values 
of using the results for gi/Fi [42]. The unpolarized structure function Fi is 
obtained from a recent parametrization of F2{x, Q"^) from NMC [43] and a fit 
to the data on Q^), the ratio of longitudinal to transverse photoabsorp- 
tion cross sections from SLAG [44]. The weight of each point is determined by 
the statistical error. The best fit coefficients are listed in Table 1. The total 
of the fits are 146 and 148 for 168 points in MS and AB schemes, respectively. 

The statistical errors on extracted parton densities Aq{x, Q^), Aq{x, Q^), and 
AG{x, Q"^) were calculated by adding in quadrature statistical contributions 
from experimental points. The weight of each point was obtained by vary- 
ing the point within its statistical error and calculating the change in the 
parton density [45]. The systematic error is usually dominated by the normal- 
ization errors (target and beam polarizations, dilution factors, etc.). Thus the 
systematic errors are to a large extent correlated point to point within one 
experiment^. We therefore assumed 100% correlated systematic errors for any 
given experiment and added systematic contributions within one experiment 
linearly. Systematic errors for each experiment were then added quadratically 
to obtain systematic uncertainties on parton densities. 

The biggest source of theoretical uncertainty comes from the uncertainty on 
the value of as- We estimate it by repeating the fits[^ with as^M^) varied 
in the range allowed by the unpolarized fixed target DIS experiments [23] 
as{Mz) = 0.108 — 0.116. The quality of the fits deteriorated significantly 
when the values as high as as{M'^) = 0.120 were used. The scale uncertainty 
is included in the error on as- We also vary current quark masses in the range 
rric = 1 — 2 GeV and mh = 4 — 5 GeV which affects the running of a^. The 
effect of SU(3)flavor breaking is estimated by varying the parameter from 1 
to 0. These factors are found to have a small influence on the results. To test 
the sensitivity to the shape of the initial distributions and the value of the 

^ We have checked that the data are consistent with this assumption. 

^ This includes both proton and deuteron data taken in a single experiment, such 

as EMS and SMC. 

^ We also relax the positivity constraints. 



8 



starting scale Qq, wc repeat the fit with initial unpolarized distributions taken 
from Ref. [1] at Qq = 1 GeV^ and find the results consistent with values given 
in Tables 1 and 2 within quoted statistical uncertainties. Possible higher twist 
effects are neglected since they are expected to drop with the photon-nucleon 
invariant mass squared as [46]. The cut > 4 GeV^ has been 

applied to all the data with the majority of them exceeding W'^ > 8 GeV^. 



4 Results and discussion 

Results for the structure functions of the proton and neutron Qi and at 

5 GeV^ are compared to the experimental data in Fig. 1. Despite a small num- 
ber of free parameters, the fits are excellent. In addition, at the initial scale 

= 0.34 GeV^ the low x behavior of the distributions is consistent with the 
Regge theory prediction ^ [20]. However, Regge theory in the past has 
been applied at the ~ 2 — 10 GeV^ of the experiments. This procedure 
clearly cannot be applied to the E154 neutron data for 0.014 < a; < 0.1, and 
is incompatible with the pQCD predictions [47,48] . If instead, Regge theory is 
assumed to apply at our starting to fix the values of powers 7/ to anywhere 
between and 1, the data are fit with four parameters. Two of those param- 
eters control the small contributions of gluons and antiquarks. When these 
parameters are fixed to zero, the resulting fit with only two free parameters 
still provides a reasonable description of the data everywhere except the low 
X region where it underestimates the E154 data on by about two standard 
deviations. 

The values of the first moments of parton distributions, as well as the first 
moments of structure functions at = 5 GeV^, are given in Table 2. The 
procedure of fitting structure functions to power laws at low = 0.34 GeV^ 
evolved up to the experimental results in low x behavior that can be 
integrated to yield the first moments. If instead, the data are fit to a power 
law in X at the average ~ 5 GeV^, a significantly bigger first moment of 
is obtained [45]. Hence our results for the first moments depend strongly 
on the assumptions that we make regarding the low x behavior. However, 
the simple assumptions that we made are attractive theoretically and have 
remarkable predictive power. 

The first moment of the deuteron structure function gf that we obtain is 
smaller than that of Ref. [9]. The reason is that our assumptions about the 

low X behavior of gi result in a contribution beyond the measured region of 
Jo^^gf dx ^ -0.014 as opposed to ^ +0.001 estimated in Ref. [9] assummg 
Regge behavior at = 3 GeV^. The first moment of gi is numerically less 
sensitive to how the data are extrapolated. 
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The first moments of the valence quark distributions are determined well, but 
the moments of the sea quark and gluon distributions are only qualitatively 
constrained. We note that the contribution of the experimental systematic 
errors to the errors on the first moments of the parton distributions is com- 
parable to the statistical contribution. The full error on the first moment of 
the gluon distribution AG is bigger than quoted in Ref. [17] despite the fact 
that the new data from E154 were added. The theoretical uncertainty is also 
quite large; it could potentially be reduced if the simultaneous analysis of the 
unpolarized and polarized data was performed (including as as one of the pa- 
rameters). It is interesting to note that at = 0.34 GeV^ the orbital angular 
momentum contribution (L^) = — 0.21q;3 is consistent with zero, i.e. hclicities 
of quarks and gluons account for most of the nucleon spin. The results of the 
fits in both MS and AB schemes are consistent within errors. The fits are 
significantly less stable in the AB scheme. Note that the values of the singlet 
axial charge Ago are essentially the same for fits in both schemes. 

The contributions from the valence quarks ^i^ncc _ (^x/18)Cg ® (A-uy + 
AMv) and sea quarks and gluons 5" ^^'^+8''^°"^ = (5/9)Cg® Ag + (1/9)Cg® AG 
to the neutron spin structure function at = 5 GeV^ are shown in Fig. 2. 
One can see that the sea and gluon contributions are larger than the valence 
contributions at x ~ 10~^. Although the sea contributions to Qi are relatively 
modest in the E154 data range x > 0.01, the strong x dependence ~ x^^'^ 
observed by E154 below a; = 0.1 is largely due the sea and gluon contributions. 
An observation of a negative value of gi at lower x and higher would provide 
direct evidence of a polarized sea. 

One may note an apparent ~ 2cr disagreement of Aq^ with the value extracted 
from the neutron beta-decay [23] Aq^ = g^ = 1.2601 ± 0.0025. This is due to 
the fact that the calculation is done in NLO, and the higher order corrections 
to the Bjorkcn sum rule arc not taken into account. The corrections can be as 
large as 5% [49] at ~ 5 GeV^. They would bring Ag3 in better agreement 
with the beta decay data. For consistency with the NLO approximation, we 
do not include this correction; it has no cfi'ect on the physical observable gi. 

Using the parametrization of the parton distributions, one can obtain the 
polarized structure function (Eq. (1)) and evolve the experimental data points 
to a common (Q^) using the formula: 

grixi, {Q')) = gr{x„ - Ag^\x„ Ql {Q')) (10) 

with 

^9f{x^, Ql m) = 9f{x^, Ql) - 9f\x,, {Q')) , (11) 
where g'^^^Xi., Ql) is the structure function measured at the experimental kine- 
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matics, and gf^ is the fitted value. The errors on gl^'^{xi, (Q^)) have three 
sources: 

<j'{gr{x,, (g^))) = a^g^U,. + a'{gr\ys,. + a\g,U,. , (12) 

where statistical and systematic uncertainties take into account the correla- 
tion between gf^^^Xi, ) and gf^, and the evolution uncertainty includes only 
uncorrelated experimental uncertainties as well as theoretical uncertainties 
added in quadrature. 

The data on the structure function g^ from two independent spectrometers 
used in E154 are given in Table 3. These results provide new information 
on the dependence of gi and thus supplement our previously published 
data [13]. Table 3 also lists the E154 data points evolved to (Q^) = 5 GeV^ 
using the MS parametrization. The NLO evolution is compared to the tradi- 
tional assumption of scaling of gi/F'^ in Fig. 3. The difference is only shghtly 
smaller than the precision of the present-day experiments. The effect on is 
small only if the integral is evaluated at the average value of (as is usu- 
ally done). The dependence of the ratio gi/ Fi is shown in Fig. 4. We plot 
the difference between the values of gi/Fi at a given Q"^ and = 5 GeV^ 
to which the SLAG data are evolved. For the neutron, the evolution of g'^ is 
slower than that of F^. Therefore, assuming scaling of g'^/F^, one typically 
overestimates the absolute value of 5'"(x, (Q^)) at low x (where Qf < (Q^)), 
and underestimates it at high x (where Qj > {Q^))- The two effects approx- 
imately cancel for the integral over the measured range in the case of E154. 
However, the shape of the structure function at low x affects the extrapolation 
to X = 0. The effect of the perturbative evolution is qualitatively the same for 
the proton. 

The data on g^ averaged between two spectrometers are given in Table 4. 
Integrating the data in the measured range, we obtain (at Q"^ — 5 GeV^) 

0.7 

J dx g^{x) = -0.035 ± 0.003 ± 0.005 ± 0.001 , (13) 

0.014 

where the first error is statistical, the second is systematic, and the third is due 
to the uncertainty in the evolution. This value agrees well with the number 
-0.036 ± 0.004 (stat.) ± 0.005 (syst.)[13] obtained assuming the indepen- 
dence of g^/F^. Using the MS parametrization to evaluate the contributions 
from the unmeasured low and high x regions, we determine the first moment 

F^ = -0.058 ± 0.004 (stat.) ± 0.007 (syst.) ± 0.007 (evol.) (14) 
at g2 ^ 5 GeVl 
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The behavior of the purely non-singlet combination {Qi — gi){x) is expected to 
be softer at low x than its singlet counterpart [50]. Evolving the E154 neutron 
and E143 proton [8] data to = 5 GeV^ and using the MS parametrization 
of Table 1 to determine the contributions from the unmeasured low and high 
X regions, we obtain for the Bjorken sum 

1 

rr"(5 GeV^) ^ j dx {gl - = 0.171 ± 0.005 ± 0.010 ± 0.006 , (15) 



where the first error is statistical, the second is systematic, and the third is 
due to the uncertainty in the evolution and low x extrapolation. This value 
is in good agreement with the 0[a\) [49] prediction 0.188 evaluated with 
as{M^) — 0.109, and it also agrees very well with the value in Table 2 obtained 
by direct integration of the parton densities. The result is fairly insensitive 
to the details of the low- a; extrapolation which is well constrained by the 
data. The low x behavior in the non-singlet polarized sector is also relatively 
insensitive to the higher-order corrections [51]. On the other hand, the low- 
X extrapolation of the proton and neutron integrals alone still relies on the 
assumption that the asymptotic behavior of sea quark and gluon distributions 
can be determined from the present data, and that the effects of higher-order 
resummations are small. These assumptions, and therefore the evaluation of 
the total quark helicity AS, are on potentially weaker grounds. Precise higher 
energy data on the polarized structure functions of both proton and neutron 
are required to determine this quantity. 



5 Conclusions and Outlook 

Additional high precision data from SLAG experiment E155 on the polarized 
structure functions of the proton and deuteron will be important in under- 
standing the spin structure of the nucleon. New results on the proton structure 
function Qi from SMG have recently been presented [52]. Also, the polarized 
electron-proton coUider experiments proposed at HERA [53] would be of great 
importance in unraveling the low x behavior of the spin-dependent structure 
function Qi. Furthermore, polarized fixed-target experiments at the Next Lin- 
ear Collider would determine the structure functions of both the proton and 
the neutron over a broad kinematic range [54], and thus compliment the HERA 
program. Extrapolations based on the fits in this paper suggest that gi{x) will 
be large at low x and have a significant dependence. Observing become 
negative at low x would provide direct evidence of a polarized sea. 

In conclusion, we have performed a Next-to-Leading order QGD analysis of 
the world data on polarized deep inelastic scattering. The data constrain the 
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first moments of tlie polarized valence quark distributions; the polarized gluon 
and sea quark distributions can only be qualitatively constrained. We deter- 
mine that the dependence of the ratio gi/Fi for the proton and neutron 
is sizable compared to present experimental uncertainties. We use the NLO 
pQCD evolution to determine the first moments of the spin dependent struc- 
ture functions of the proton and neutron at = 5 GeV^, and find that the 
data agrees with the Bjorken sum rule within one standard deviation. 
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Table 1 

Fitted values of the free parameters in Eq. (6) in MS and AB schemes. Also quoted 
are the statistical, systematic, and theoretical errors. 





MS 


AB 




Value 


Stat. 


Syst. 


Theory 


Value 


Stat. 


Syst. 


Theory 







99 


+0.08 
-0.08 


+0.04 
-0.05 


+0.97 
-0.11 





98 


+0.07 
-0.06 


+0.05 
-0.07 


+0.96 
-0.09 


Ad 


-0 


78 


+0.14 

-0.20 


+0.05 

-0.05 


+0.05 

-1.28 


-0 


82 


+0.06 
-0.11 


+0.07 

-0.06 


+0.31 

-1.21 


Aq 


-0 


02 


+0.03 
-0.06 


+0.01 
-0.02 


+0.01 
-0.35 


-0 


04 


+0.02 
-0.05 


+0.01 
-0.02 


+0.03 
-0.06 


Ag 


1 


6 


+1.1 
-0.9 


+0.6 
-0.6 


+0.2 
-1.3 





1 


+1.0 
-0.3 


+0.5 
-0.2 


+0.1 
-0.6 


au 





63 


+0.06 
-0.07 


+0.04 
-0.05 


+0.36 
-0.06 





55 


+0.08 
-0.06 


+0.03 
-0.04 


+0.56 
-0.05 







28 


+0.15 
-0.11 


+0.05 
-0.03 


+0.75 
-0.03 





40 


+0.20 
-0.12 


+0.07 
-0.13 


+0.53 
-0.34 


aq 





04 


+0.29 
-0.03 


+0.12 
-0.03 


+0.55 
-0.01 





00 


+0.17 
-0.00 


+0.17 
-0.00 


+0.00 
-0.00 


ac 





8 


+0.4 
-0.5 


+0.3 
-0.3 


+0.1 
-0.6 








+0.7 
-0.0 


+ 1.0 
-0.0 


+1.0 
-0.0 



Table 2 

First moments of the polarized parton distributions and structure functions of the 
proton, neutron, and deuteron in MS and AB schemes evaluated at = 5 GeV^. 
Errors are statistical, systematic, and theoretical. 









MS 






AB 






Value 


Stat. 


Syst. 


Theory 


Value 


Stat. 


Syst. 


Theory 







69 


+0.03 
-0.02 


+0.05 
-0.04 


+0.14 
-0.01 


0.74 


+0.03 
-0.02 


+0.03 
-0.03 


+0.07 
-0.01 


Adv 


-0 


40 


+0.03 
-0.04 


+0.03 
-0.03 


+0.07 
-0.00 


-0.33 


+0.03 
-0.04 


+0.03 
-0.05 


+0.01 
-0.03 


AQ 


-0 


02 


+0.01 
-0.02 


+0.01 
-0.01 


+0.00 
-0.03 


-0.03 


+0.02 
-0.02 


+0.01 
-0.01 


+0.01 
-0.01 


AG 


1 


8 


+0.6 
-0.7 


+0.4 
-0.5 


+0.1 
-0.6 


0.4 


+1.0 
-0.7 


+0.9 
-0.6 


+1.1 
-0.1 


Aq3 


1 


09 


+0.03 
-0.02 


+0.05 
-0.05 


+0.06 
-0.01 


1.07 


+0.03 
-0.02 


+0.06 
-0.06 


+0.10 
-0.01 


Aqs 





30 


+0.06 
-0.05 


+0.05 
-0.05 


+0.23 
-0.01 


0.42 


+0.05 
-0.08 


+0.06 
-0.06 


+0.03 
-0.01 


AS 





20 


+0.05 
-0.06 


+0.04 
-0.05 


+0.01 
-0.01 


0.25 


+0.07 
-0.07 


+0.05 
-0.05 


+0.05 
-0.02 


Ago 





20 


+0.05 

-0.06 


+0.04 

-0.05 


+0.01 

-0.01 


0.21 


+0.05 

-0.06 


+0.06 
-0.07 


+0.05 
-0.02 


F? 


0.112 


+0.006 
-0.006 


+0.008 
-0.008 


+0.009 
-0.001 


0.114 


+0.005 
-0.006 


+0.010 
-0.011 


+0.001 
-0.003 


^ 1 


-0.056 


+0.006 
-0.007 


+0.005 
-0.006 


+0.002 
-0.001 


-0.051 


+0.005 
-0.006 


+0.006 
-0.007 


+0.001 
-0.012 


ri 


0.026 


+0.005 
-0.006 


+0.005 
-0.006 


+0.005 
-0.001 


0.029 


+0.004 
-0.005 


+0.007 
-0.008 


+0.001 
-0.007 




0.168 


+0.005 
-0.004 


+0.008 
-0.007 


+0.007 
-0.001 


0.165 


+0.004 
-0.004 


+0.009 
-0.009 


+0.013 
-0.001 
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Table 3 

E154 results on gf at the of the measurement for each spectrometer. Also shown 
are results for (/" evolved to (Q^) = 5 GeV^ according to Eq. (10). Errors were 
propagated as described in the text. 





Ql 




5f(xi,5GeV2) 




GeV2 


lb Stat. lb syst. 


± Stat. ± syst. ± evol. 


2.75° spectrometer 


0.017 


1.2 


-0.351 ±0.115 ±0.109 


-0.421 ± 0.115 ± 0.113 ± 0.016 


0.024 


1.6 


-0.374 ± 0.071 ± 0.064 


-0.409 ± 0.071 ± 0.066 ± 0.007 


0.035 


2.0 


-0.289 ± 0.061 ± 0.038 


-0.304 ± 0.061 ± 0.039 ± 0.005 


0.049 


2.6 


-0.212 ±0.041 ±0.022 


-0.215 ± 0.041 ± 0.023 ± 0.004 


0.078 


3.3 


-0.119 ±0.031 ±0.013 


-0.117 ± 0.031 ± 0.013 ± 0.002 


0.123 


4.1 


-0.075 ± 0.030 ± 0.010 


-0.073 ± 0.030 ± 0.010 ± 0.001 


0.173 


4.6 


-0.070 ± 0.033 ± 0.010 


-0.069 ± 0.033 ± 0.010 ± 0.001 


0.241 


5.1 


-0.053 ± 0.028 ± 0.008 


-0.053 ± 0.028 ± 0.008 ± 0.000 


0.340 


5.5 


0.002 ± 0.036 ± 0.004 


0.001 ± 0.036 ± 0.004 ± 0.000 


0.423 


5.8 


0.027 ± 0.059 ± 0.007 


0.027 ± 0.059 ± 0.007 ± 0.000 


5.5° Spectrometer 


0.084 


5.5 


-0.152 ±0.029 ±0.019 


-0.153 ± 0.029 ± 0.019 ± 0.001 


0.123 


7.2 


-0.117 ±0.017 ±0.013 


-0.121 ± 0.017 ± 0.013 ± 0.002 


0.172 


8.9 


-0.059 ± 0.016 ± 0.009 


-0.066 ± 0.016 ± 0.009 ± 0.003 


0.242 


10.7 


-0.040 ± 0.012 ± 0.006 


-0.047 ± 0.012 ± 0.006 ± 0.003 


0.342 


12.6 


-0.019 ± 0.012 ± 0.005 


-0.024 ± 0.012 ± 0.005 ± 0.001 


0.442 


13.8 


-0.009 ± 0.012 ± 0.003 


-0.011 ± 0.012 ± 0.003 ± 0.001 


0.564 


15.0 


0.003 ± 0.008 ± 0.001 


0.003 ± 0.008 ± 0.001 ± 0.000 
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Table 4 

Combined results on at the of the measurement. Also shown are results for 
gi evolved to (Q^) = 5 GeV^ according to Eq. (10). 



Xi 








5f(xi,5GeV2) 




GeV2 


lb Stat. 


lb syst. 


± Stat. ± syst. ± evol. 


0.017 


1.2 


-0.351 ± 


115 ±0.109 


-0.421 ± 0.115 ± 0.113 ± 0.016 


0.024 


1.6 


-0.374 ± 


071 ± 0.064 


-0.409 ± 0.071 ± 0.066 ± 0.007 


0.035 


2.0 


-0.289 ± 


061 ± 0.038 


-0.304 ± 0.061 ± 0.039 ± 0.005 


0.049 


2.6 


-0.204 ± 


040 lb 0.022 


-0.207 ± 0.040 ± 0.023 ± 0.004 


0.081 


4.5 


-0.137 ±0 


021 ±0.016 


-0.136 ± 0.021 ± 0.016 ± 0.002 


0.123 


6.6 


-0.108 ±0 


015 ±0.012 


-0.111 ± 0.015 ± 0.012 ± 0.002 


0.173 


8.2 


-0.061 ± 


014 ± 0.009 


-0.067 ± 0.014 ± 0.009 ± 0.003 


0.242 


9.8 


-0.042 ± 


Oil ± 0.007 


-0.048 ± 0.011 ± 0.007 ± 0.003 


0.342 


11.7 


-0.017 ±0 


Oil ± 0.005 


-0.021 ± 0.011 ± 0.005 ± 0.001 


0.441 


13.3 


-0.007 ± 


Oil ± 0.002 


-0.009 ± 0.011 ± 0.002 ± 0.001 


0.564 


15.0 


0.003 ± 


008 ± 0.001 


0.003 ± 0.008 ± 0.001 ± 0.000 



Q = 5 Ge\' 
• SLACE143 
o CERNSMC 
— OurNLOMt 




-0.0« 



Q = 5 GeV 
• SLACE154 
o CERNSMC 
— OurNLOFit 




10 



10 



10 



10 



Fig. 1. The structure functions xc/l and xgf at = 5 QeV^. The E143, SMC, 
and E154 data have been evolved to = 5 GeV^ using a procedure described in 
the text. The result of the MS fit is shown by the solid line and the hatched area 
represents the total error of the fit. 



18 




10 10 10 

X 



Fig. 2. The contributions to the strueture function of the neutron from the 
valence quarks [(1/18)C5 ® (Auy + 4A(iv)] (sohd hne) and from the sea quarks and 
gluons [(5/9)C5(8)AQ + (1/9)Cg<8)AG] (dashed hne). The shaded and hatched areas 
represent the total uncertainties on each quantity. 
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0.01 



n 



0.005 



-0.005 



-0.01 



-0.015 



-0.02 



-0.025 



10 



-2 



= 5 GeV^ E154 data 



• NLO evolution 
O gj/Fj scaling 



Systematic uncertainty 



Evolution uncertainty 



10 
X 



-1 




Fig. 3. The structure function xqi evolved to = 5 GeV^ using our 
parametrization and using the assumption that Qi/Fi is independent of Q^. 
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Proton 



Neutron 




5 10 15 20 

(GeV^) 



5 10 15 20 

(GeV^) 



Fig. 4. Evolution of the ratios gi/Fi for proton (left) and neutron (right). Plotted is 
the difference f^l^' *3^) ~ |^(^) 5 GeV^). The MS fit is shown by the solid line and 
the hatched area represents the total (experimental and theoretical) uncertainty of 
the fit. 
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